Non-Perturbative Functional Renormalization Group for Random Field Models and 
Related Disordered Systems. I: Effective Average Action Formalism 



o 
o 

(N 
o 

Q 

o 



Si 

O 



in 



c 

o 
o 



> 
o 
in 
in 

cn 

(N 

o 



X 



Gilles Tarjufl and Matthieu TissieiQ 

LPTMC, CNRS-UMR 7600, Universite Pierre et Marie Curie, 
boite 121, 4 PI- Jussieu, 75252 Paris cedex 05, France 
(Dated: February 2, 2008) 

We have developed a nonperturbative functional renormalization group approach for random field 
models and related disordered systems for which, due to the existence of many metastable states, 
conventional perturbation theory often fails. The approach combines an exact renormalization group 
equation for the effective average action with a nonperturbative approximation scheme based on a 
description of the probability distribution of the renormalized disorder through its cumulants. For 
the random field O(N) model, the minimal truncation within this scheme is shown to reproduce the 
known perturbative results in the appropriate limits, near the upper and lower critical dimensions 
and at large number N of components, while providing a unified nonperturbative description of the 
full (N, d) plane, where d is the spatial dimension. 
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I. INTRODUCTION 

The effect of quenched disorder on the long-distance 
physics of many-body systems largely remains an unset- 
tled question despite decades of intensive research. On- 
going controversies persist for instance on the equilibrium 
and out-of-equilibrium behavior of spin glasses and sys- 
tems coupled to a random field*^ 2 - Even though progress 
has been made, it has so far proven difficult to construct 
a proper renormalization group (RG) approach provid- 
ing a description of ordering transitions and criticality in 
these systems. A technical reason for this unsatisfactory 
situation is that quenched disorder makes the system in- 
trinsincally inhomogeneous and that one should in prin- 
ciple follow the renormalization of the whole probability 
distribution of the disorder. A physical reason is that the 
presence of disorder and of the resulting spatial inhomo- 
geneity lead, for at least some range of the control pa- 
rameters, to multiple "metastable states" . (At this point 
we use the term "metastable state" in a loose acceptance 
to describe configurations that minimize some energy or 
free-energy, action or effective action in field-theoretical 
terminology, but differ from the true ground state.) How 
such metastable states evolve upon coarse-graining under 
RG then represents the central issue: at large length- 
scale, their influence could vanish, leaving only benign 
signatures in the thermodynamics, or else it could mod- 
ify the critical behavior of the system, the nature of its 
phases, and, often in an even more spectacular way, the 
relaxation and out-of-equilibrium dynamical properties. 

A well known example of the kind of puzzles associ- 
ated with quenched disorder and metastable states is the 
failure of the so-called "dimensional reduction" property 
in the random field Ising model (RFIM).— d&H Standard 
perturbation theory predicts to all orders that the critical 
behavior of the RFIM in dimension d is the same as that 
of the pure Ising model, i.e., in the absence of random 
field, in two dimensions less, d — 2. The property has 
been shown in a compact and elegant manner by Parisi 



and Sourlas 7 - by means of a supersymmetric formalism. 
However, dimensional reduction predicts a lower critical 
dimension for ferromagnetism in the RFIM of rf/ c = 3, in 
contradiction with rigorous results!^ The dimensional 
reduction property must therefore break down in low 
enough dimension. The supersymmetric approach gives 
a hint at the origin of the breakdown, which appears to 
be related, yet in a somewhat obscure way, to the pres- 
ence of multiple metastable states^ (in this case, local 
minima of the Hamiltonian) . 

Over the years, and on top of numerous computer sim- 
ulations and scarce exact analytical results, theoretical 
approaches have been devised to cope with disordered 
systems characterized by multiple metastable states, such 
as spin glasses and random field models. 1 To list the main 
ones, we mention: 

(i) phenomenological approaches such as the heuris- 
tic domain-wall argument o 11 ' 12 and the "droplet" 
description ; 13 ' 14 i 15 in which one directly focuses on rare 
excitations and the associated low-energy metastable 
states; 

(ii) mean-field theories, combined with the replica for- 
malism in order to handle the average over disorder; for 
models with spin-glass ordering, the potentially dramatic 
effect of the metastable states is captured through a spon- 
taneous breaking of the replica symmetry 1 6 1 1 7 

(iii) specific RG techniques for low-dimensional (d = 
1,2) systems, as for instance the Coulomb gas RG ap- 
proach for two-dimensional disordered XY model s 18 ' 19 
or real space RG for strongly disordered one-dimensional 
systems ; 20 ! 21 ' 22 

(iv) the perturbative functional RG for energy- 
dominated disordered models considered in the vicinity 
of a critical dimension at which the fundamental fields 
are dimensionless ; 23 ' 24 '^^'^ 7 -'^ one must then follow the 
flow of a whole function, an appropriate renormalized cu- 
mulant of the disorder. As shown first by Fisher— for 
an elastic manifold pinned by a random potential, the 
long-distance physics is controlled by a zero-temperature 
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fixed point at which the renormalized cumulant is a non- 
analytic function of the fields, with the nonanalyticity 
encoding the effect of the many metastable states at zero 
temperature. 

All these approaches, however, are either question- 
able or not easily generalizable: on the one hand, the 
phenomenological approaches lack rigorous foundations 
and the relevance of mean-field descriptions to finite- 
dimensional systems is, to say the least, far from garan- 
teed; on the other hand, the perturbative functional 
RG becomes extremely complex, and soon untractable 
in practice for random field systems when going beyond 
one-loop calculations j 29 i 30 ' 31 moreover, it does not allow 
one to study the RFIM (as for specific RG techniques, 
they are not extendable by construction). 

The purpose of the present work, described here and 
in a companion paper— is therefore to propose a gen- 
eral theoretical framework that leads to a consistent de- 
scription of the equilibrium behavior of the random field 
models and related disordered systems. To achieve this, 
we rely on a version of Wilson's continuous RG via mo- 
mentum shell integration, 33 Under various terminologies, 
"Exact RG", "Functional RG", and "Nonperturbative 
RG" , it has been developed in the past 15 years to become 
a powerful method for investigating both universal and 
nonuniversal properties in Statistical Physics and Quan- 
tum Field Theory, 3 ^^! 3 -^ The approach is "exact" 
in the sense that the RG flow associated with the pro- 
gressive account of the field fluctuations over larger and 
larger lengthscales is described through an exact func- 
tional differential equation. It is "functional" because 
through the exact equation, one follows the flow of an 
infinite hierarchy of functions of the fields in place of sim- 
ply coupling constants. It is "nonperturbative" (beyond 
the mere tautology that an exact description automati- 
cally includes all perturbative as well as nonperturbative 
effects) because it lends itself to efficient approximation 
schemes that are able to capture genuine nonperturbative 
phenomenal to name a few, in the case of the standard 
O(N) scalar model, (numerically) tractable approxima- 
tions describe the Kosterlitz-Thouless transition of the 
XY model in d = 2, known to be associated with the 
binding/unbinding of topological defects (vortices), as 
well as the convexity property of the thermodynamic po- 
tential in case of spontaneous symmetry breaking, recov- 
ered in other treatments through nonperturbative config- 
urations like instantons. 

To study the problem at hand, we combine the 
ideas of the perturbative functional RG for disor- 
dered systems with the general formalism of the ex- 
act /functional/nonperturbative RG. In the following, we 
shall denote our approach nonperturbative functional RG 
(NP-FRG). It provides a framework to study both per- 
turbative and nonperturbative effects in any spatial di- 
mension d and for any number of components of the 
fundamental fields, N. We exclude from the scope 
of the present series of articles relaxation and out-of- 
equilibrium dynamic phenomena, as well as spin glass 



ordering. We also postpone to a forthcoming publica- 
tion the development of the NP-FRG in a superfield 
formalism able to directly address the failure of super- 
symmetry in connection with that of dimensional reduc- 
tion. Short versions of the present work have appeared 
in Refs. 



The present paper is organized as follows. 

In section II we present the models and the formalism. 
We first introduce the models and discuss their physical 
relevance and the main open questions. From the corre- 
sponding replica field theories, we then derive the exact 
RG equation for the effective average action, which is the 
generating functional of the one-particle irreducible cor- 
relation functions at the running scale. We next relate 
the replica formalism, in which the replica symmetry is 
explicitly broken through the application of sources, to 
the cumulants of the renormalized disorder. We close the 
section by writing down the exact RG flow equations for 
these cumulants. 

In section III, we introduce a systematic nonpertur- 
bative approximation scheme. After first discussing the 
symmetries of the problem and the way to implement 
them in the effective average action formalism, we intro- 
duce the nonperturbative truncation scheme of the exact 
RG equation: it relies on (i) an expansion in cumulants 
of the disorder and (ii) a well tested approximation of the 
nonperturbative RG, the "derivative expansion" , which 
uses the fact that the relevant physics is dominated by 
long wavelength modes to perform an expansion in the 
number of spatial derivatives of the fundamental fields. 
Finally, we detail the minimal truncation that we use 
in our numerical investigation of the random field O(N) 
model (RFO(AT)M). 

In section IV, we specialize the formalism to the study 
of the RFO(7V)M. We introduce the scaling dimensions 
suitable to a search for the putative zero-temperature 
fixed point controlling the ordering transition. We first 
consider the case of the RFIM and then extend our de- 
scription to the RFO(A r )M. With the help of these di- 
mensions, the RG flow equations are then cast in a scaled 
form. We also briefly comment on possible application 
to other disordered systems. 

We next discuss in section V an important property of 
the truncations previously described: because of the one- 
loop structure of the exact flow equations and of the ap- 
propriate choice of the approximations, one recovers the 
perturbative results both near the upper critical dimen- 
sion, d uc = 6, and in the N — > oo limit of the RFO(iV)M. 
Even more interestingly, we also show that our minimal 
truncation near the lower critical dimension for ferromag- 
netism of the RFO(N > 1)M, d tc = 4, reduces to the 
perturbative functional RG result (at one loop) obtained 
from the nonlinear sigma version of the models To the 
least, the truncated NP-FRG thus provides a nonper- 
turbative interpolation in the whole (N, d) plane of the 
known perturbative results near d = 4, d = 6, as well as 
N -> oo. 

Finally, the presentation and discussion of the results 
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obtained for the RFO(N)M within the present NP-FRG 
approach will be described in the companion paper ^ 



are both taken for simplicity (see also discussion below) 
with gaussian distributions characterized by zero means 
and variances given by 



II. MODELS AND FORMALISM 
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A. Models 

We focus on the equilibrium, long-distance behavior 
of a class of disordered models in which iV-component 
classical variables with O(N) symmetric interactions are 
coupled to a random field. Depending on whether the 
coupling is linear or bilinear, the models belong to the 
"random field" (RF) or the "random anisotropy" (RA) 
subclasses. Such models with N = 1, 2, or 3 are relevant 
to describe a variety of systems encountered in condensed 
matter physics or physical chemistry. To name a few, one 
can mention dilute antiferromagnets in a uniform mag- 
netic field^ critical fluids and binary mixtures in aero- 
gels (both systems being modelled by the N = 1 RF 
Ising model) ; 41 i 42 ' 43 vortex phases in disordered type-II 
superconductors (described in terms of an elastic glass 
model whose simplest version is the N = 2 RF XY 
model) ) 44 i 45 ' 46 amorphous magnets, such as alloys of rare- 
earth compoundS ) 47 i 48 and nematic liquid crystals in dis- 
ordered porous media (described by A^ = 2 or N = 3 RA 
models) 4^ 

Other related models can be described as well within 
the same formalism, but will only be alluded to: the "ran- 
dom elastic" model describing an elastic system, such as 
an interface or a vortex lattice, pinned by the presence of 
impurities; the "random temperature" model associated 
with impurity-generated bond or site dilution in a ferro- 
magnetic Ising model. For reasons that will become clear 
further down in this section, we exclude from the present 
study spin glass ordering and we rather concentrate on 
ferromagnetic ordering (in which the O(N) symmetry is 
spontaneously broken) or "quasi-ordering" (phases with 
quasi-long range order) . 

Our starting point is the field theoretical (coarse- 
grained) description of the systems in terms of an N- 
component scalar field xi x ) m a (^-dimensional space and 
an effective Hamiltonian, or bare action, 



Jx *• Z u=l 



A. 



V 



N . 

where J = J d d x and the superscript \x spans the N 
components of the field; h(x) is a random magnetic field 
and t(x) a second-rank random anisotropy tensor, which 



T^{x)TP°{y) = — {5^5^ + S ( x - V). ( 3 ) 

where the overbar generically denotes the average over 
quenched disorder. Higher-order random anisotropics 
could be included as well. They will indeed be gener- 
ated along the RG flow. However, for symmetry reasons, 
when starting with only a second-rank, or more gener- 
ally an even-rank, random anisotropy, only even-rank 
anisotropics are generated: this corresponds to what is 
called the random anisotropy (RA) model. The model 
with a nonzero A, for which anisotropics of both odd and 
even ranks are generated under RG flow, is the random 
field (RF) model. 

The equilibrium properties of the model are obtained 
from the average over disorder of the logarithm of the 
partition function, 

Z[J] = J P X cxp (-S[x;h,r}+ J J(x)-x(x) \ , 

(4) 

where J(x) is a source linearly coupled to the fundamen- 
tal field and a (ultra-violet) momentum cutoff A, associ- 
ated with an inverse microscopic lengthscale such as a 
lattice spacing, is implicitly considered in the functional 
integration over the field. With this definition however, 
the partition function and the corresponding thermody- 
namic potential W[J] = lnZ[J] are still functionals of 
the random fields: W[J] = W[J; h, r]. As is well known 
from the theory of systems with quenched disorder, the 
thermodynamics is given by the average over disorder of 
the "free energy", i.e., 



W[J] = lnZ[J]. 



(5) 



Full information on the system, in particular an access 
to the correlation (Green) functions of the field, requires 
knowledge of the higher moments of W[«7], viewed as 
a random functional. 74 As will be discussed more thor- 
oughly further below, such information can be conve- 
niently extracted by using the replica formalism whose 
starting point is the replacement of InZ by the limit of 
(Z n — l)/n when n, the number of replicas of the original 
system, goes to zero. Quite differently from the standard 
but controversial use of this replica trick, in which the 
analytic continuation for n < 1 opens the possibility of a 
spontaneous breaking of the replica symmetry^ we will 
consider an a priori more benign procedure in which the 
symmetry between replicas is explicitly broken by the in- 
troduction of external sources acting on each replica in- 
dependently. This procedure will allow us to generate the 
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cumulant expansion of the disorder-dependent functional 
W[J}. 

Within the replica formalism, the original problem 
is replaced by one with n replica fields {Xa(x)}, a — 
1,2, ■•• , n, and the "replicated action", obtained after 
explicitly performing the average over the disorder in the 
partition function, reads: 

Sn[{Xa}}= I l\j2[\dXa{x)\ 2 +T\ X a{ X )\ 2 + 
Jx L Z a=l 

1 " 

^(\Xa(x)\ 2 ) 2 ]-- J2 [ A Xa(x) ■ \b(x) 



6=1 
2 



12 vi/v»v n , j 2 

+ A 2 ( Xa (x) ■ Xb(x)) 
with the corresponding partition function 

Zn [{Jo}] = / f[ V X a exp ( - S n [{ X a}} 
J a=l ^ 

n p 

+ J2 J *(x)-Xa(x] 

1 J X 



(6) 



(7) 



where the linear sources J a (x), a = 1,2,-- - ,n, act on 
each replica separately. Associated to this partition func- 
tion is the generating functional of the connected Green 
functions, W„[{J a }] = lnZ„[{J a }], and the effective ac- 
tion, r„[{0 a }], defined through a Legendre transform: 



n - 

r„[{0a}] = -W n [{ Ja}} + £ / Ja(x) • (f> a (x), (8) 



the fields {<fi a } and the sources {J a } being related by 

6W n [{J a }} 



r a (x) = (xS(sb)) 



6J£(x) 



(9a) 



where (X) represents the average of X with the weight 
given m Eq. Q, and 



(5^(05) 



(9b) 



The effective action is the generating functional of the 
one-particle irreducible (1 — PI) correlation functions or 
proper vertices. 

The formalism we are about to describe also applies to 
extensions of the replicated action of Eq. © that can be 
cast in the form 



Sn [{Xa}} =/{E [\\dXa{x)f + U A ( X a(x))} 
Jx *• o=l Z 

1 " 

- 9 £ V A { Xa (x), X b(x)) + -- 



a, 6=1 



(10) 



where the subscript A recalls that the various terms are 
at their bare value, defined at the microscopic scale A, 
and the dots indicate possible functions involving higher 
numbers of replicas. The functions U\, V\, ■ ■ ■ satisfy the 
O(N) symmetry as well as the S n permutational symme- 
try between replicas. Eq. |T]) is obviously a special case 
of the above expression, and higher-order anisotropics 
are included in a 2-replica term which is only function of 
Xa( x ) ■ Xb{x). RF and RA O(N) models with nongaus- 
sian distributions of the random fields and anisotropics 
are described by terms involving higher number of repli- 
cas. (Note that the RA O(N) model is defined as such 
for A > 1; the Ising case, N = 1, corresponds to an- 
other model, the random temperature one introduced 
hereafter.) 

Other disordered systems are also described by the 
form of the replicated action in Eq. (fT0|) . For instance, 
the random temperature model corresponds to Eq. (|10|) 
with U\ and Va functions of the fields only through the 
O(N) invariants p a = \\Xa\ 2 Pb = ^\Xb\ 2 - In the RF, 
RA, and random temperature models, the 1-replica part 
of the bare action simply describes n copies of the stan- 
dard ferromagnetic O(N) model without disorder. 

The random elastic model is also a special case of 
Eq. (fT0| . However, contrary to the models just discussed, 
the 1-replica potential U\ is absent (or reduced to a 
purely quadratic term), so that there is no mechanism 
triggering a paramagnetic-ferromagnetic phase transi- 
tion. The 2-replica potential Va, which is the second 
cumulant of a random pinning potential, is now func- 
tion of only the difference between the two replica fields, 
Xa{x) — Xb(x). As a result, the model has an additional 
symmetry, the statistical tilt symmetry— which garan- 
tees that the 1-replica part of the action, including the 
kinetic term, is not renormalized: the effective action has 
thus the same 1-replica part as the bare one. (Note that, 
as shown in Ref. [301 ] and in the companion paper ^ the 
random elastic model, albeit with an underlying period- 
icity, also emerges as a low-disorder approximation of the 
RF and RA XY (N = 2) models.) 



Exact RG equation for the effective average 
action 



The exact RG in the effective average action 
formalis m 34 ' 36 i 51 relates the bare action, here Eq. (jlOp . 
to the full effective action, Eq. (jSJ), through a progressive 
inclusion of fluctuations of longer and longer wavelength. 
To do so, one introduces an infrared regulator, charac- 
terized by a scale fc, which, in the functional integration 
leading to the partition function, suppresses the contri- 
bution of the low-energy modes with momentum |q| < k 
while including the high-energy modes with \q\ > k. Af- 
ter Legendre transformation, this defines an "effective 
average action" at the running scale fc, T/j, which contin- 
uously interpolates between the microscopic scale k — A, 
at which Tj~=a reduces to the bare action, and the macro- 
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scopic one, k = 0, at which T k= o equals the full effective 
action. 

More precisely in the present context, a "mass-like" 
quadratic term is added to the bare action, Eq. (JTUJ) , 

, n N r 

AS fc [{* Q }] = o E E / K: ab il 2 )x^-q)xl(q), 

a, 6=1 n,v=l J 1 

(11) 

where f q = J d d q/(2n) d ; R k v ab {q 2 ) denotes infrared cut- 
off functions which, in order to enforce that the additional 
term satisfies the same 0(N) and S n symmetries as the 
bare action (see above), must take the following form: 

K V M) = (Mq 2 )s ab + Rk(q 2 )) <V (12) 

The cutoff functions R k (q 2 ) and Rk(q 2 ) are chosen 
such as to realize the decoupling of the low- and high- 
momentum modes at the scale k: for this, they must 
decrease sufficiently fast for large momentum |q| ^> k 
and go to a constant value (a "mass" ) for small momen- 
tum \q\ <C k. The presence of an off-diagonal component 
Rk(q 2 ) is somewhat unusual and will be discussed later 
on. The cutoff functions must also satisfy the two con- 
straints that (i) they go to zero when k — ► 0, so that one 
indeed recovers the full effective action with all modes ac- 
counted for, and (ii) R k {q 2 ) diverges while Rk{q 2 ) stays 
finite when k — > A, so that the effective average action 
does reduce to the bare action. (In what follows we are 
only concerned with the long-distance behavior of the 
models and do not pay attention to microscopic details; 
we thus let A go to oo in the cutoff functions.) Differ- 
ent choices have been proposed and tested in the recent 
literature. Standard choices for R k (q 2 ) are of the form 



R k (q 2 ) = Z k q 2 r(q 2 /k 2 ) 



(13) 



with Z k a field renormalization constant yet to be spec- 



ified and r{y) 



y 



l (l - y)Q(l - y)^- where 9 is the 



Heaviside function, or r(y) = (e v — l)" 1 ^! 

From the partition function Z k [{J a}] obtained from 
the bare action supplemented with the fc-dependent reg- 
ulator, Eq. (fTTj) . one defines the generating functional 
of the Green functions Wfe[{J }] = lnZfe[{J Q }] and, 
through a Legendre transform, one has access to the ef- 
fective average action at the running scale k, T k : 

n „ 

r k [{4> a }] + W k {{J a }} =J2 / JaN-0aN-A<S fe [{0 a }] 

(14) 

where the fields {<fi a } and the sources {J a } are related 
by the (fc-dependent) expression 

# (*) = <x£(«)> = ^ 

OJa(X) 

The Legendre transform is slightly modified by the ad- 
dition of the last in Eq. (fH)) , which ensures that the 



effective average action T k does reduce to the bare ac- 
tion at the microscopic scale, with no contribution from 
the infrared regulator. This addition does not change 
the behavior in the k — ► limit since the regulator goes 
identically to zero. Physically, and to use the language of 
magnetic systems, the effective average action is a coarse- 
grained Gibbs free energy. It is the generating functional 
of the 1 — PI correlation functions from which one can 
derive all Green functions of the modified system at the 
scale k. Note that here and in the following we omit the 
subscript n associated to the number of replicas in order 
to simplify the notations. 

The evolution of the effective average action with the 
infrared cutoff k is governed by an exact flow equation, 



Tr\d k R k {q z 



,(2) 



Ri 



(16) 

where the trace involves a sum over both replica indices 
and A^-vector components; R k (q 2 ) is defined in Eq. (fT12|) 

(2) 

and T k is the tensor formed by the second functional 
derivatives of T k with respect to the fields <f)% (q) : 



\ / ab 



S 2 T, 



(17) 



The above RG flow equation is a complicated func- 
tional integro-differential equation that cannot be solved 
exactly in general; but, due to its one-loop structure and 
its reasonably transparent physical content, it provides 
a convenient starting point for nonperturbative approxi- 
mation schemes. 

At this point, it is quite clear to see why we have ex- 
cluded spin glass ordering from our considerations. The 
quadratic form of the infrared regulator in Eq. (1111) sup- 
presses the fluctuations of the low-momentum modes of 
the fundamental fields Xa- Spin glass ordering on the 
other hand involves fluctuations of composite fields, asso- 
ciated e.g. to the "overlap" between different replicas. 16 
Proper RG treatment of such fluctuations implies to in- 
troduce a "mass-like" regulator for composite fields, i.e., 
in the simplest case a functional that is quartic in the fun- 
damental fields instead of the quadratic term used here. 
We do not consider this case in the present work. 



C. Explicit replica symmetry breaking and 
cumulants of the renormalized disorder 

Among the technical difficulties encountered when 
making use of the exact RG equation, Eq. (HHJ), there 
is one which is specific to disordered systems and to the 
present replica formalism: one must invert the matrix 

(2) 

F y k ab + R ky ab for arbitrary replica fields (since all replicas 
are different due to the independently applied sources). 
Before delving into this problem, it is worth giving some 
physical insight into the meaning of the explicit replica 
symmetry breaking used here. 
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As discussed in section II- A, after full account of 
the fluctuations, the bare disorder is renormalized to 
a full random ("free energy") functional W[J], which, 
to make its dependence on the bare quenched disorder 
explicit, we now denote W[J;/i]. This random object 
can be characterized by the infinite set of its cumulants, 
W l [ Jr] , W 2 [ Ji , J 2 ] , W 3 [ J u J 2 , J 3 ] , • • • , with 



W 1 [J 1 ]=W[J 1 ;h] 



(18) 



W 2 [J U J 2 ] = W [J i ; h] W[ J 2 ; h] - W [Ji ; h] W [ J 2 ; h] , 

(19) 

etc... The first cumulant W\ gives access to the thermo- 
dynamics of the system and the higher-order cumulants 
describe the distribution of the renormalized disorder (we 
define, as in the bare action, a disorder with zero mean). 
Note that by construction the cumulants are invariant 
under permutations of their arguments. 

The cumulants can be generated from an average in- 
volving copies, or "replicas" , of the original disordered 
system, as follows: 



a=l 



exp(J2 w l J *> = exp (W[{ J a }}) 

(n 1 n 

]TWMJa]+ g E W 2l J a> J b] (20) 
a=l a, 6=1 

n 

E W 3 [J a ,J b ,J c } + --- 



a,b,c—l 

where the n copies have the same bare disorder but are 
coupled to different external sources. To fully character- 
ize the random functional W[J; h], it is indeed important 
to describe its cumulants for generic arguments, i.e., for 
different sources. (Be aware that the subscripts 1,2,... 
used to denote the cumulants of W should not be con- 
fused with the subscript n denoting the number of repli- 
cas in section II-A and omitted since: here for instance, 
Wi denotes the 1-replica component, corresponding to 
the first cumulant, whereas with the previous notation 
W n —i is given by the sum of all cumulants with all there 
arguments equal.) 

A convenient trick to extract the cumulants with their 
full functional dependence is to let the number of replicas 
be arbitrary and to view the expansion in the right-hand 
side of Eq. (|20[) as an expansion in increasing number of 
"free" , or unconstrained, sums over replicas of the func- 
tional W[{ J a}] defined below Eq. ([7]). The term of order 
p in the expansion is a sum over p replica indices of a 
functional depending exactly on p replica sources, this 
functional being precisely equal here to the pth cumu- 
lant of W[J;/i]. This procedure, which rests on an ex- 
plicit breaking of the replica symmetry and an analytic 
continuation to arbitrary numbers of replicas (including 
the limit n — > previously introduced), is a priori differ- 
ent from the standard use of replicas, in which all sources 



are equal, and it avoids the delicate handling of a spon- 
taneous replica symmetry breaking 16 ' 1 7 It has been 
used in a similar context by Le Doussal and Wiesei 53 i 54 
The practical implementation of the expansion in free 
replica sums will be detailed in the next subsection. 

In our present NP-FRG approach however, the central 
object is the effective action F, not W. The expansion of 
r[{</> a }] in increasing number of free replica sums reads 



n n 

r[{«M] = E r ^«] - o E r 2[^,0 6 ] 



.6=1 



^| E F3[<i> a ,<i>b,<j>c] H 



(21) 



a,6,c=l 



where for later convenience we have introduced a minus 
sign for all even terms of the expansion. r[{0 a }] and 
14^ [{J a }] are related by a Legendre transform, so if one 
also expand the sources J a [{4>f}} (where we have denoted 
{4>f} the n replica fields to avoid confusion in the indices) 
in increasing number of free replica sums, one can relate 
the terms of the expansion of the effective action to the 
cumulants of the random functional W[J;fe]. The rela- 
tion is straighforward for the first terms, but gets more 
involved as the order increases. 

More precisely, Ti [(f)] is the Legendre transform of 
Wi[J], namely, 



T 1 [d)} = -W 1 {J]+ f J(x)-4>(x), 



with 



5Wi[J] 



and the second-order terms is given by 

r»[<h.,<h] = w 2 [j[<h],j[4>2]], 



(22) 



(23) 



(24) 



where J[(f>] is the nonrandom source defined via the in- 
verse of the Legendre transform relation in Eq. (|22p . i.e., 
J"[0](ac) = Sr 1 [4>]/Scj^(x). (Note that J(x) introduced 
here differs from the source J a (x) introduced in equa- 
tion (|9b[) : through the Legendre relations, the latter de- 
pends on all the fields {<p a } while the former depends on 
only one replica field.) The above expression motivates 
our choice of signs for the terms of the expansion in free 
replica sums of r[{0 a }], Eq. (f2"0")) : ^[cfrijfa] is directly 
the second cumulant of W[J; h] (with the proper choice 

Ofj[0]). 

For the higher-order terms, one finds 



r 3 [01,02 ) 03] = -W 3 [J[cf> 1 ],J[cf> 2 ],J[cl> 3 }} + 

K 2 Vw>i]])~ J 

\ / xv 



(25) 



x W$[J[<h], J[<h]] +perm(123) 
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etc., where perm(123) denotes the two additional terms 
obtained by circular permutations of the fields 4>i , cf> 2 , <j>% 
and where we have used the following short-hand nota- 
tion: 



w\ p) 

1,031 . 



W. 



(pi) 



2,x 1 ...x p ,y 1 ...y, 



.1* 



6Ji(xi)...6Ji(x p )' 



5P+iW 2 [J u J 2 ] 



(26) 



5Ji{xi)...8Ji(x p )8J 2 {yi)...8J 2 (y q ) ' 



(27) 

etc. Note that for clarity the O(N) indices have been 
omitted in the above expressions. 

We point out that T p [<pi, <p p ] for p > 3 cannot be 
directly taken as the pth cumulant of a physically accessi- 
ble random functional, in particular not of the disorder- 
dependent Legendre transform of W[J; h] (although it 
can certainly be expressed in terms of such cumulants 
of order equal or lower than p) . In the following and by 
abuse of language, we will nonetheless generically call the 
T p s "cumulants of the renormalized disorder" (which is 
true for p = 2). 

In complement to the above picture and more specif- 
ically for random field systems, it is also interesting to 
introduce a renormalized random field (or random force) 
h[cf>](x) defined as the derivative of a random free-energy 
functional, 



D. Exact RG equations for the renormalized 
disorder cumulants 



The reasoning developed in the previous subsection can 
be applied to the effective average action Tk and its ex- 
pansion in free replica sums. As a results, Eqs. (|18H29p 
can be extended to any running scale k. Yet, to make 
the expansion in free replica sums an operational proce- 
dure, one needs be able to perform systematic algebraic 
manipulations, as for instance the inversion of the ma- 
trix appearing in the right-hand side of the exact RG 
equation, Eq. (fTo) . We detail here the method for matri- 
ces depending on two replica indices, but functionals of 
the n replica fields. Extension to higher-order tensors is 
presented in Ref. [5^ |. 

A generic matrix A o &[{0/}], where we have again de- 
noted {4>f} the n replica fields to avoid confusion in the 
indices, can be decomposed as 



(30) 



In the above expression, it is understood that the second 
term A a b no longer contains any Kronecker symbol. Each 
component can now be expanded in increasing number of 
free replica sums, 



X[{0/}]=l [ ° ] [0a]+^l [1I [0a|0c] 



(31) 



h[<j>r( x ) 



b»(x) 



(w[J[<f>};h]-W[J[<f>};h]) , (28) 4»&[{0/}] = AM[<f> a ,<j> b ] + Y j A^[<j> a ,<j> b \<j> c 



and whose first moment is equal to zero by construc- 
tion. It is easy to derive that its pth cumulant (p > 2) 
is given by the derivative with respect to <pi,...,<p p of 
Wp[J[0i], J[0 P ]], which can then be related to deriva- 
tives of T2,r3, for instance, 



h[cp 1 ](x)h[ ( f )2 }(y) = r 2 n J y [ ( jy 1 ,4>- 



(29) 



(32) 

where the superscripts in square brackets denote the or- 
der in the expansion (and should not be confused with su- 
perscripts in parentheses indicating partial derivatives). 

(2) 

As an illustration, the expansion of the matrix 1^ 
defined in Eq. (|17[) reads, in terms of the expansion of 
effective average action itself, 



where have used a short-hand notation similar to that of 
Eqs. (|26l27p and omitted the TV- vector indices for sim- 
plicity. Terms of order 3 and higher are again given by 
more complicated expressions. 

We close this discussion by noticing that in the sim- 
pler case of the random manifold model, Ti and W\ being 
trivial and unrenormalized due to the statistical tilt sym- 
metry (see above), J[4>] has a simple explicit expression. 
For instance, if the bare action has a quadratic 1-rcplica 
term, Y\ [<p] is equal to this quadratic functional and J[4>] 
is a known linear functional of </>, which further simpli- 
fies when considering uniform fields. This allows one to 
devise ways to directly measure the second cumulant of 
the renormalized disorderi 55 ' 56 Nothing similar occurs in 
random field and random anisotropy models: the thermo- 
dynamics of such systems being highly nontrivial (with 
a phase transition and a critical point) , the expression of 
J[cj>] is involved and a priori unknown. 



ri 2) [{0/}] Q = r^[0 a ]- 



I 



E 1 ^ 



[4>a,<f>c 



r£ 2) [{*/}]«* = -r& ; 



_ _-r( n ) 



E 



,(110) 
fe.3 



[4> a , 4> h , 4> c 



(33) 



(34) 



where the permutational symmetry of the arguments of 
the I\ p 's has been used. 

Algebraic manipulations on such matrices can be per- 
formed by term-by-term identification of the orders of 
the expansions. For instance, the inverse B = A^ 1 of 
the matrix A can also be put in the form of Eq. ((30|) 
and its components, B a and B a b, expanded in number of 
free replica sums. The term-by-term identification of the 
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condition A ■ B = 1 leads to a unique expression of the 
various orders, and B p \ of the expansion of B in 
terms of the A^ 's and A^ 's with q < p. The algebra be- 
comes rapidly tedious, but the first few terms are easily 
derived: 



m[ ( f> 1 ]=A^[ ( f> 1 }-' 



(35) 



BM&ufc] = -BM[<t>i]AM[<f>i,<p 2 ]m[<p 2 ] (36) 



£M[0i|0 2 ] = -BM^IW [0i |</> 2 ]£ [0] [0i] (37) 



£?W [01 , 02 10 3 ] = -B [01 [0l] | AW [0! , 2 10 3 ] 
-i[°][01,0 3 ]S[ O ][03]l[°][0 3 ,0 2 ] 

-lW[0i|0 3 ] J B[°][0i]i[ o ][0i,0 2 ] 

- JP [0i , 2 ]Bl°] [0 2 ] II 1 ! [0 2 |0 3 ] IbM [0 2 



(38) 



etc. 

We can apply the above procedure to the exact RG 
equation for the effective average action. For conve- 
nience, we introduce the modified propagator at the scale 
k, 



Pk[{<f>f}} 



,(2) 



+ Rk 



with 



Pfc,oh[{0/}] = Pk,a[{<t>f}] S ab + Pk,ab[{4>f}\, 



(39) 



(40) 



where P k ,a and P k , a b are still tensors with respect to 
momenta and vector component indices. Eq. (|16p then 
leads to an infinite hierarchy of flow equations for the 
cumulants of the renormalized disorder, 

d k T kA [0i] = ~ / ld k (R k {q 2 ) +Rk(q 2 ))trpM_ q [0i] 



d k R k (q*)trpM_ q [0i,0 



(41) 



d k T k>2 [0i,0 2 ] = — J [d k (R k (q 2 ) + R k (q 2 )) 

trP ~ki- g [<h\<h] + d k R k {q 2 )trP [ ^ q _ q [01,0x102] 
+ d k R k {q 2 )trP [ ^_ q [0 1 ,0 2 ]+p e rm(12)|, 

(42) 

and so on, where tr indicates a trace over TV-vector com- 
ponents and perm(12) denotes the expression obtained 
by permuting 0i and 2 . (Some care is needed in the 



term by term identification in order to properly sym- 
metrize the expressions and satisfy the permutational 
property of the various arguments of the cumulants.) 

Expressing the higher-order terms Pu and P^ with 
p > 1 only by means of P| ' and the derivatives of the 
Tjt^'s and introducing the short-hand notation d k to in- 
dicate a derivative acting only on the cutoff functions, 
i.e., d k = d k R k S/SR k + d k R k 5/5R k , Eq. (JUJ can be 
rewritten 

d k Y K2 [0i, 2 ] = ±d k Trlpf ] [0i] (if 2 0) [0 l5 2 ] 



,(110) 



[0!, 01, 2 ]) + Pf 1 [01, 0!] if 2 0) [0!, 2 ] 



~pf [0i, 2 ] (rj$ [0i, 2 ] - R k l) +perm(12) 



(43) 



and similarly for higher-order cumulants , where l qq , — 
{2-K) d 8{q + q')5^ u and the trace Tr is now over both mo- 
menta and and N- vector components; the modified prop- 
agators Pj^ and P] ' are explicitly given by 



Pi o] [0i] = (r^[0i] + i? fe i 



(2) | 



(44) 



if ] [0 l5 2 ] = Pf [0l](r^ J [01, 2 ] - R k l)Pr [02] 

(45) 

This provides a hierachy of exact RG equations for 
the cumulants of the renormalized disorder (including 
the first one which leads to a description of the ther- 
modynamics). One should note that (i) the cumulants 
are functional of the fields and contain full information 
on the complete set of 1 — PI correlation functions and 
(ii) the flow equations are coupled, the (p + l)th cumu- 
lant appearing in the right-hand side of the equation for 
the pth cumulant. As such these RG equations remain 
untractable and their resolution requires approximations. 



_ BP] 



(ii) [ 



>[0][ 



III. NONPERTURBATIVE APPROXIMATION 
SCHEME 

A. Symmetries in the effective average action 
formalism 

When writing the RG flow for the effective average ac- 
tion and when devising an approximation scheme to solve 
it, one should as far as possible make sure that the sym- 
metries of the theory are not explicitly violated at any 
scale. Such a requirement is easily implemented as far 
as elementary symmetries, such as invariance by trans- 
lation and rotation in Euclidean space, O(N) symmetry, 
and S n replica permutational symmetry, are concerned: 
the infrared regulator AS k added to the bare action must 
be chosen such that it is invariant under the appropriate 
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transformations, which is indeed garanteed by the expres- 
sions in Eqs. (|11I12[) . The exact effective average action 
at any scale k then also possesses the symmetries of the 
bare action, and one just had to be careful that the trun- 
cations do not explicitly break the symmetries, which is 
easily implemented^ 

A similar treatment can be applied to most additional 
symmetries of the disordered systems under considera- 
tion. For instance, the "statistical tilt symmetry" of 
the random manifold model is easily extended to a k- 
dependent statistical tilt symmetry with any regulator 
of the form given in Eqs. (| 1 111 2|) . which implies that 
the 1-replica part (first cumulant) of the effective av- 
erage action is unrenormalized along the flow. Simi- 
larly, the additional inversion symmetries of the random 
anisotropy (x a • Xb — > ~Xa ■ Xb) and the random tem- 
perature (xo, {Xb} b ^ a -> ~Xa, {Xb\bfa) models are read- 
ily accounted for with the choice Rk = 0. Truncation 
schemes naturally follow. 

Taking into account the underlying supersymmetry 
that characterizes the random field model for a gaussian 
distribution of the random field 7 is much more involved. 
First, because one knows that the supersymmetry, which 
goes with the dimensional reduction property, must be 
broken in low enough dimension (at least, in d — 3), so 
that, even if the RG flow is started with an initial con- 
dition obeying supersymmetry, a mechanism should be 
provided to describe a spontaneous breaking of the su- 
persymmetry. Secondly, the supersymmetry shows up in 
a superfield formalism built with auxilliary fermionic and 
bosonic fields, but it is far from transparent in the present 
framework based on the fundamental fields. (This is true 
already at the level of the initial condition of the RG 
flow.) We shall therefore defer the proper resolution of 
this problem to a forthcoming publication.— Note that 
an underlying supersymmetry is also present in the ran- 
dom manifold model, where it also leads to the d — > d— 2 
dimensional reduction. However, the pure model with 
no disorder is merely a free field theory, and this is easily 
accounted for^ 



B. Truncation schemes 



truncation scheme is provided by the so-called "deriva- 
tive expansion" i 34 ' 36 It consists in expanding the effec- 
tive average action in increasing number of derivatives 
of the field(s) and retaining only a limited number of 
terms. The lowest order is the "local potential approxi- 
mation" (LPA) 5 4 in which one only considers the flow of 
the effective average potential, i.e., the effective average 
action for a uniform field configuration. The field is not 
renormalized and the associated anomalous dimension is 
equal to zero. Field renormalization, which is important 
in the present problem where one expects the anomalous 
dimension to be quite sizeable in low dimensions (e.g., 
numerical estimates give n ~ 0.5 for the RFIM in d = 3), 
requires to go beyond the LPA and to consider the first 
order of the derivative expansion. Previous studies on a 
variety of systems, including the pure O(N) model, have 
shown that the system's behavior is quantitatively very 
well described at this level of approximationi 36 i 37 ' 60 i 61 
Higher-order terms improve the accuracy ; 62 ! 63 but they 
rapidly become untractable except in simple models. 

For the disordered systems considered here, one more 
step is needed. We have seen in section II- C that an 
expansion in number of free replica sums can be used 
to generate the cumulants of the renormalized disorder. 
Keeping only a limited number of terms in the expansion 
therefore leads to a systematic truncation scheme. To 
describe both the thermodynamics and the renormalized 
probability distribution of the disorder, one must con- 
sider at least the first two cumulants, or equivalently, the 
second order in the expansion in free replica sums. 

Finally, on top of the two previous approximations, it 
may be useful, and numerically more tractable, to ex- 
pand the functions appearing in the truncated effective 
average action in powers of the field considered around 
a given (uniform) configuration. This configuration can 
be taken either as zero everywhere or as a nontrivial con- 
figuration that minimizes the effective average potential 
(here, more precisely, its 1-replica component that gives 
access to the thermodynamics). Again, the accuracy and 
convergence properties of such field expansions have been 
widely tested for many different models. In the present 
case, and for reasons that will become clear later on, field 
expansions should be used with great caution. 



We have already stressed that solving the exact RG 
equation for the effective average action requires approx- 
imations. The general framework has proven quite ver- 
satile for devising efficient and numerically tractable ap- 
proximations which are able to describe both universal 
and nonuniversal properties in any spatial dimension and 
to capture genuine nonperturbative phenomena (see In- 
troduction). Such approximations generally amount to 
truncating the functional form of the effective average 
action, which results in a self-consistent flow that pre- 
serves the fundamental structure of the theory (as the 
symmetries, see above). 

If one is interested in the long-distance physics of a sys- 
tem and in observables at small momenta, a systematic 



C. Minimal truncation 

Given the general scheme presented above, the choice 
of a minimal nonperturbative trucation is guided by a 
combination of factors: experience gained from studies 
on other models, constraints associated with the symme- 
tries of the full theory, intuition or previous knowledge 
concerning the physics of the problem at hand, require- 
ment of being able to recover as much as possible exact 
and perturbative results in the appropriate limits, and of 
course, a practical limitation coming with the numerical 
capability to actually solve the set of RG flow equations. 

As we have already alluded to, a description of the 
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long-distance physics of random field models and related 
disordered systems at least requires to keep the first two 
cumulants of the disorder, i.e., the first two terms, T^i 
and I\2, of the expansion of the effective average action 
in free replica sums. Because of the anticipated nonnegli- 
gible value of the anomalous dimension of the field rj, one 
must also include in the description at least the first or- 
der of the derivative expansion of the first cumulant 
The resulting truncated functional form of the effective 
average action then reads 



r fc [{*„}]= / IE 

Jx •> a=l 



Uk(Pa(x)) + ~Z k ( Pa (x))\dcp a (x)\ 2 



-Y k (p a (x))(dp a (x)) 2 



1 " 1 
- 2 E V k (<l> a (x),cp b (x))\, 

a,b=l ' 



(46) 



where, as before, p a (x) = ^\<p a (x)\ 2 . In the above ex- 
pressions, Uk((pi) = Uk{pi) is the effective average poten- 
tial, which is equal to the 1-replica component T^i eval- 
uated for a uniform field and will hereafter be simply de- 
noted the 1-replica potential; V k {<p\, <t>2) = V k ( P i, p2, 4>i ' 
<fi 2 ) is the 2-replica potential and is equal to the 2-replica 
component IV2 evaluated for a uniform field configura- 
tion. Physically, U k ((pi) is a coarse-grained Gibbs free 
energy and V k {(pi,(p 2 ) is the second cumulant of the 
renormalized disorder evaluated for uniform fields (see 
Eqs. (HI [HI)). The two terms Z k (pi) and Y k {p x ) cor- 
respond to field renormalization functions for the Gold- 
stone and massive modes, respectively. 

We note in passing that the fact that only the first two 
cumulants of the disorder have been kept in the trunca- 
tion does not imply that the probability distribution of 
the renormalized disorder is actually taken as gaussian. 
Indeed, as will be discussed in the companion paper^ 
the probability is not gaussian in general. The trunca- 
tion means that we have neglected the contribution com- 
ing from the third cumulant in the RG flow of the second 
cumulant and have therefore decoupled the hierarchy of 
flow equations for the cumulants. 

Being interested in the description of the models in the 
full (N, d) diagram, we will have recourse to further ap- 
proximations that make the numerical resolution of the 
flow equations easier. More specifically, we consider the 
lowest-order term of the field expansion of the field renor- 
malization functions around a nontrivial configuration, 
Pm.k — 5- 1 0m,fc | 2 3 which minimizes the 1-replica potential 
U k (p): Y k = and Z k (p) = Z m ^ k , with Z m ^ k = Z k (p m . k ) 
and U k ( Pm> k) = 0. Physically, (p m ^ k is the magnetiza- 
tion (order parameter) at the scale k. (If cp m ,k^o = 0, 
the system is in an O(N) symmetric phase whereas if 
4>m,k— >o 7^ 0, the system is in the phase with broken 
symmetry.) Z m ^ k is chosen as the field renormalization 
in the cutoff function R k (q 2 ) (see Eq.JT3j)). 

Finally, we simplify the resulting RG flow equations by 
setting the off-diagonal cutoff function to zero, R k = 0. 
As will be shown, this choice leads in general to an 



explicit breaking of dimensional reduction (despite the 
fact that the infrared regulators vanish identically when 
k — ► 0). In the following paper— we shall discuss the way 
to nonetheless make sense out of the results, the distinc- 
tion between spurious and real breaking of dimensional 
reduction being easily characterized. A complete reso- 
lution of this issue will be provided when extending the 
NP-FRG approach to the superfield formalism^! 

With the above approximations which we shall refer 
to as the minimal trucation, the self-consistent NP-FRG 
equations can be derived from Eqs. (|4"lTf4"5)) . The flows 
of the 1- and 2-replica potentials read 



1 



dkU k (pi) = ~ J d k R k (q 2 )tr^Pl° ] (q 2 - lPl ) 

-Pl \ q 2 -,p 1 )V k ^(<l> 1 ,c Pl )P k 0] (q^ Pl : 



d k V k {cp l ,cp 2 ) = -- J d k R k (q 2 )tr\Pf\ q 2 - Pl ) 
V^\< Pu 4> 2 ) + P k 0] (q 2 -,p 2 )V k {11 \< Pu <p 2 ) 2 + 



(47) 



vf o) (0 lj 3 )Pr(8 a ;Pi)^"'(^i.0x) 
V^ 1 \c Pu cl> 1 )Pl \ q 2 ;p 1 )V^\^c P2 ) 



(H), 



+ perm(12) >, 



(48) 

where the trace is over the N- vector components and, due 
to the O(N) symmetry, V k (<p l7 <f> 2 ) = V k ( Pl ,p 2 ,z) with 

z = 4>i-cp 2 /y/'ipip 2 ; the (modified) propagator P k ^ (q 2 ; p) 
is given by 



(1 - <W) 



Z m , k q 2 +R k (q 2 ) + U' k {p) 

SjA 

Z m ,kq 2 + Rk(q 2 ) + U' k (p) + 2pU' k '{p) 



(49) 



where p, = 1 is chosen to be the direction of the field <fi 
and therefore corresponds to the massive mode while the 
(N — 1) remaining components represent the Goldstonc 
modes. 

The flow of the field renormalization constant 
Z mk is obtained from the prescription Z k (p) = 

d q 2T k 2 \(q 2 ]pY Ltl \ q 2 =0 with p chosen as a Goldstone mode 
{ji 7^ 1)2& and from the condition U' k {p m _ k ) = 0. It can 
be explicitly written as 



d k Z m , k — d a 2d k 



(50) 



q=0 
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where a line denotes the Goldstone propagator and dots 
represent vertices obtained from derivatives of either the 
1-replica potential (single dots) or the 2-replica potential 
(dots linked by a dashed line); for instance, 




(21) (21) 

represents the three-point vertex 1^ = V£ . We did 
not include the graphs containing 4-point vertices be- 
cause in the truncation considered here, they do not con- 
tribute to the flow of Z m ^- From the above flow equa- 
tion, Eq. (|50p . one extracts a running anomalous expo- 
nent, 

i] k = -kd k Z m ^ k . (51) 



The initial conditions for the RG flow equations are 
obtained from the bare action, Eq. (fT0)l . The RG 
flow equations form a closed set of coupled nonlinear 
integro-differential equations for two functions, U k (pi) 
and Vk(pi, P2, z), and a constant, Z rn ±. The numerical 
task of solving these equations is still arduous and when 
needed for reducing the difficulty of the computations, 
we will also consider truncated expansions of the 1- and 
2-replica potentials in some or all of their field arguments 
(see below). 

The present approach represents a nonperturbative 
but of course approximate RG description. Already at 
the minimal truncation discussed above, one includes 
all operators previously suggested to be important for 
capturing the long-distance behavior of the present dis- 
ordered models, namely operators involving 1- and 2- 
replica terms. As will be shown further below, it also 
reduces to the leading results of perturbative RG analy- 
ses near the upper critical d imension, d uc — 6, near the 
lower critical dimension for ferromagnetism when N > 1, 
d = 4, and when the number of components N becomes 
infinite. One of its main advantages is that it provides a 
unified framework to describe models in any spatial di- 
mension d and for any number N of field components. As 
such, it garantees a consistent interpolation of all known 
results in the whole (N, d) plane, in addition to allow- 
ing the study of genuine nonperturbative phenomena. If 
more accuracy is needed, the truncation scheme proposed 
in III-B gives a systematic means to refine the descrip- 
tion, by including e.g. the third cumulant or a more 
detailed account of the momentum dependence of the 
1 — PI vertices. 



IV. RANDOM FIELD MODEL 

A. Scaling dimensions near a zero-temperature 
fixed point 

For the RFIM, it has been propose d 64 i 65 , and con- 
vincingly supported by numerical and experimental 
results ^ 40 ' 66 , that the fixed point controlling the crit- 
ical behavior associated with the transition between a 
high-temperature - or large-disorder strength - disor- 
dered (paramagnetic) phase and a low-temperature - or 
small-disorder strength - ordered (ferromagnetic) phase 
is at zero temperature (see Figure [1]). The existence of 

A A 




FIG. 1: Schematic phase diagram of the RFIM in the disorder 
strength A - temperature T plane above the lower critical di- 
mension di c = 2 (temperature can be introduced at the bare 
level through the Boltzmann weight). At low disorder and low 
temperature, the system is ferromagnetic, and it is paramag- 
netic otherwise. The arrows describe how the renormalized 
parameters evolve under the RG flow at long distance, and 
/ and RF denote the critical fixed points of the pure and 
random-field Ising models, respectively. 

such a zero-temperature fixed point around which tem- 
perature is dangerously irrelevant leads to a somewhat 
anomalous scaling at the critical pointi 64 ' 65 The two in- 
dependent critical exponents characterizing the scaling 
behavior of the pure Ising model should a priori be sup- 
plemented by an additional exponent 9 describing the 
vanishing of the (renormalized) temperature as the fixed 
point is approached. This exponent 9 leads to a modi- 
fication of the so-called hyperscaling relation, which be- 
comes 2 — a = [d — 9)v where the critical exponents a 
and v have their usual meaning, and to a new scaling 
of the correlation functions. In particular, the so-called 
"connected" and "disconnected" components of the pair 
correlation function (or 2-point Green function) behave 
at the critical point as: 

G c (q) = (x(-q)x(q)) - (x(-q))(x(q)) ~ (52) 



In the following, we more specifically focus on the ran- 
dom field O(N) model. 



G d (q) = (x(-q))(x(q)) ~ (53) 
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where rj is the usual anomalous dimension of the field and 
fj is related to the temperature exponent 9 through 



fj = 2 



(54) 



Above the upper critical dimension d uc — 6, the expo- 
nents take their classical, mean-field values, rj = 0, a = 
0, v = 1/2, and 9 = 2, leading to fj = 0. The dimensional 
reduction property leads to a constant shift of dimen- 
sion, d — > d — 2, i.e., to 9 = 2 and fj = r\, all exponents 
being in addition given by those of the pure model in di- 
mension d— 2. Whether the scaling behavior around the 
critical point is described by 3 independent exponents, or 
only 2, has been a long-time issue, with suggestions that 
an additional relation applies, 9 = 2 — rj or equivalently 
fj = 2r/£2- We shall address and answer this question in 
the following paper i^ 2 . 

To search for a zero-temperature fixed point, it is con- 
venient to introduce a renormalized temperature. Ac- 
tually, one could add an explicit temperature T in the 
Landau-Ginzburg- Wilson description of the model con- 
sidered here: multiplying the argument of the exponen- 
tial in the partition function, Eq. by a factor T" 1 
to make the correspondence with the Boltzmann factor 
of Statistical Physics leads to a bare replicated action in 
Eqs. ([6]) and (| 10[) in which the 1-replica part, including 
the kinetic term, is multiplied by a factor T , the 2- 
replica part by T~ 2 , etc. Generally speaking, one can use 
this temperature T as a book keeping device to sort the 
orders in the expansions in number of free replica sums. 
As a result for instance, the modified propagator P^ [</>i] 

is proportional to T whereas Pj^ [<pi , 2 ] is independent 
of T. One can use this book keeping trick to devise ways 
to define a renormalized temperature at running scale k, 
Tk, which reduces to the "bare" temperature T at the 
microscopic scale k = A. To this end, we first define the 
renormalized disorder strength at scale k, A m ,fc, as 



4>m,k, 02 — 0m, fc) ■ 



(55) 



where as before 4> m ,k is a field configuration that mini- 
mizes the (1-replica) potential U k (<f>), and Afc(0i,02) is 
the second cumulant of the renormalized effective ran- 
dom field defined as in Eq. (f2"9l . namely, 

A fc (0i,0 2 ) = y fe (11) (0i,0 2 ). (56) 

In the present truncation, the second cumulant is only 
considered for homogeneous field configurations and rj! 1 ^ 
reduces to V k with the same notations for partial 

derivatives as in Eqs. (23 123) (e.g., V fc (U) (0i, 2 ) = 
dcj) 1 d^ 2 Vk (0i,02)). At the microscopic scale A, A m ,k re- 
duces to Aa/T 2 where A a is the bare variance of the 
random field and the factor T~ 2 comes for reasons just 
explained above. 

A running temperature can now be defined by 



T k = 



Z m ,k (fc/A) 2 
(A m u/A A f 



(57) 



One checks that since Z m ^ = T _1 (see Eq. (fTOj) and 
discussion above), Tk indeed reduces to T when k = A. 
An associated running exponent is obtained from 



= fcd fc lnT fe . 



(58) 



By using the definition of r/k, one may alternatively in- 
troduce a running exponent fjk = 2 — 9k + r/k, which 
converges to the critical exponent fj defined in Eqs. (|531 
l54| if the relevant fixed point is reached, and compute it 
from the equation 



r/k - 2r)k = kdk&r. 



(59) 



On top of the usual scaling dimensions, Uk, 14 ~ k d and 
~ C^m k& rf ~ 2 ) 1/ ' 2 ) one can use the running temperature 
to define dimensionless quantities (denoted by lower-case 
letters) suitable for looking for a zero-temperature fixed 
point: 



Zm.k Tk 



1/2 



Uk (0l) = 7f u k (<Pl) 
+ k 



k d 

V k (01,02) = 7f^V k (ipi,<P2) 
1 k 



±k 



n. 



-S k {ipi,ip 2 ) , 



(60a) 
(60b) 
(60c) 
(60d) 



with 5k (<pi, 1P2) — v k ('■Pi, ¥2)- Note that with the def- 
initions of A m , fc and T k , 5 m , k = 4 (<Pm,k,<Pm,k) is con- 
stant along the RG flow and equal to its initial value 
Aa/A 2 (in practice, and since we are not interested here 
in making a precise connection to the microscopic scale, 
we will set S m ^k = !)■ 



Scaled form of the exact RG equations for the 
RFIM 



With the use of the above defined dimensionless renor- 
malized quantities, the flow equations can be expressed 
in a scaled form. Specifically, one can recast Eqs. (|4"T|) 
and JUl for N = 1 in the form 



d t Uk(ip) = -(d - 2 + rjk - r]k)u k (<p) 
+ ^(d- 4 + fj k )ipu' k ((p) 



(61) 



+ 2v d {l[ d) K(<p)) 5 k (<p, <p) + m^iuliip))} 
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d t v k ((pi,<p 2 ) = ~(d - 4 + 2% - 2r} k )v k (<pi, <p 2 )+ 
^(d- 4 + f] k )(ip 1 d ipi +^^2)^(^1,^2) 
-2v d {l{%Uvi),<(<P2))8 k (<Pu<p 2 ) 2 

+ 4 d) «(^))4(vi^i)4 20) (^^2) ( 62 ) 



+ 4 d) K(^))4 02) (^,^)]} 



where d t is a derivative with respect to t = ln(fc/A), a 
prime denotes a derivative with respect to the field (when 
only one argument is present), v^ 1 = 2 d+1 TT d ^ 2 T(d/2), 

and we recall that 5 k ((pi, y> 2 ) = (ifi, if 2 ); In (to) 

and zi, 1 ^ra 2 (wi, w 2 ) are the "dimensionless threshold func- 
tions" defined from the infrared cutoff function, Eq. (fT3)) , 



as 



.36,51 



,(d), \ l r 1 x \ f°° j d/2Vkr(y) + 2yr'(y) 

(63) 

1 f°° 

l&nAvuv*) = - 2 / *// /2 (%r(») + 2yr'(y)) 



(p(y) + wi) ni (p(y) + w 2 )™ 2 \p(y) + ^i p(y) + w 2 



"2 



(64) 



with p(y) = j/(l + r(y)) and y = q 2 /k 2 . The proper- 
ties of these threshold functions, whose detailed behav- 
ior depends on the choice of the infrared cut-off function 
r(y), have been extensively discussed . 36 ' 51 They decay 
rapidly when w ^> 1, which, since u'{.(ip) = U' k ((f))/(Z k k 2 ) 
is the square of a renormalized mass, ensures that only 
modes with mass smaller than k contribute to the flow 
in Eqs. (fBTj) and (f6"2")l . As an illustration, the use of the 
so-called "optimized" cut-off function r(y) — y~ x (l — 
y)0(l — y), 5 - leads to explicit expressions, namely, 



1 



d + 2j (1 + wi)"i(l + w 2 )™2 



1 + W\ 1 + w 2 



(65) 



The threshold functions essentially encode the nonper- 
turbative effects beyond the standard one-loop approxi- 
mation. Note that, although not shown in the notation, 
the threshold functions explicitly depend on the scale k 
via the running exponent rj k . 

The above flow equations for u k (<Pi) and v k (ipi,ip 2 ) 
are supplemented by equations for r\ k and rj k , i-e., for 
Z m>k and T k or A m> fe. (Note that the equation for fjk 
is actually redundant as it is a consequence of the other 
equations; it is nonetheless convenient to introduce and 



use it.) The flow equation for Z m>k follows from Eq. ([50l) 
and one finds: 



,11/2 



ilk = ^{4m^« lfe ,< life )u mfc - 

2mW« fe ,< fe )< fe C, fe ( 66 ) 

where we have used the short-hand notation 5' k (ip) = 

d v S k ((p,(p) — 8fr(tp,ip) + 5^ (<p,<p) and the subscript 
"m, fe" indicates that the functions are evaluated for fields 
equal to <p mt k'i we have also introduced the additional 
(dimensionless) threshold function 



m ( nln 2 (wi,W 2 ) 



2 Jo 



1 



(p(y) + w 1 ) n *(p(y)+w 2 ) na 



dyy d/2 (l + r(y)+yr'{y)) 



(l + r(y) +yr'{y)) 



(Vkr{y) + 2yr'(y)) 



"2 



My) + w i p(y) + u, 2 

2Vk(r(y) + yr'(y)) - 4y(2r'(y) + yr" (y)) \, 



(67) 



whose properties are discussed in Ref. [36j|51(. For 
instance, with the "optimized" regulator introduced 
above— one finds that 



\tn 2 (Wl,W 2 ) (1 + Wi) „ l(1 + W2)ll2 - 



1 



(68) 



Finally the flow equation for A m , k (or equivalently the 
flow of the constraint S m> k — 1 discussed below Eq. (f6"0"l) 1 
leads to the following equation: 

2%-% = 2vJl{ d \u'^ k )uZ 2 k - J<X,k 



+ 4 \ u m.k)(Kn,k + n^m,k 



-s 



X'2 

+ '1 \ u m ,k)— 1 k 

U m,k 



J l [ u m,k)\^°m,k Tl °m,k + o m ' k ' 
z u m.k z 



1 



l 2 \ u m,k) a m.k°m.k 



m.k h 



(69) 



where, as before, 5' k ((p) = d v 5 k (ip,ip) and similarly for 
5fc(ip), and we have introduced 

Sfc(yi) = Am (d Vl - d V2 ) 2 (S k ((pi,ip 2 ) ~ S k (ipi,ipi)) 2 

(70) 



and 



EfcO^i) = lim (d Vl - d V2 ) 2 8 k (<fi,ip 2 ) . (71) 



All other notations are as before. 

Before extending the results to the KFO(N)M, we 
point out important features of the above equations. 
First, we have kept terms proportional to T k but, pro- 
vided one reaches a fixed point with an exponent 9 = 
9k— >o > where temperature is thus irrelevant, those 
terms are subdominant in the scaling region k — > 0. In 
particular, the fixed point is attained by following the 
flow with an initial temperature T equal to zero. 

Secondly, "anomalous" terms, and Tj.S m ^, ap- 

pear in the expression of 2r\ k —f\ k ■ As can be inferred from 
Eqs. (I70p and (|7T|) . Yi m ,k can only differ from zero, and 
S m .fc become infinite, when a non-analyticity (a "cusp") 
in {if>i — ip 2 ) appears in the (dimensionless) renormalized 
disorder function S k ((pi,ip 2 ) when if2 — > <p\ (and both 
go to (p m ,k)- If 8k(<pi,tp2) is analytic, no signature of 
such anomalous behavior is found. (We have implicitly 
assumed that no stronger nonanalyticity appears, which 
means that a fixed point can be reached and that the 
theory is renormalizable; this has to be checked in ac- 
tual computations.) We shall come back in more de- 
tail to these two important aspects of the NP-FRG ap- 
proach in the following paper ^2, Finally, one may notice 
that because of the Z 2 = 0(1) symmetry, the potential 
Uk is an even function of ip and because of the addi- 
tional permutation symmetry, v k (ipi, P2) = Vk(V2) fi) — 

Vk(-(pi,-<P2) = Vk(-<P2, ~<Pl)- 



C. Generalization to the RFO(iV)M 

The preceding treatment can be extended to the 
RFO(7V)M. The variable p = ^\4>\ 2 is written in terms 
of a dimensionless variable, p = k d ~ 2 T k ~ 1 Z^ l 1 k p , where 
the tilde will be dropped in the following when no confu- 
sion is possible between dimensionless and dimensionful 
quantities. The variable z = <pi ■ <fi 2 / (2^/pipli) is already 
dimensionless. 

For the 1-replica second-order tensors (in TV-vector 
components) evaluated for a uniform field configuration, 



e.g., for if 1 (q 2 ■ X ) or for A k {<p U <j>i) = ^ ll; (0 1 ,0 1 ), 
the O(N) symmetry reduces the number of terms to a 
"longitudinal" component (corresponding to the massive 
mode, see Eq. (|49jl ) and N — 1 identical "transverse" 
components (corresponding to the Goldstone modes, see 
Eq. ([4"9")0 . We therefore introduce 

"SjTO 9 ' P' z=l ) = <W [ 5 ^k,h{p) + (1 - S»i)$k, t(p)} , 

(72) 

with 



,(H)/ 



8k,L{p) = 2pd Pl d P2 v{pi,p 2 ,z = l)\ Pl , 



Sk,r(p) = —d z v(p, p, z)\ z= x, 
2p 



(73) 



(74) 



and we define the longitudinal, uik,L{p), and transverse, 
w k,r{p), masses as 

w kiL (p)=u' k {p) + 2pu' k '(p), (75) 



W k ,T(p) = u' k (p), 



14 

(76) 



where a prime now denotes a derivative with respect to 
P- 

The renormalized disorder strength at the running 
scale k can be characterized, e.g., through the transverse 
component, A ki T(p, P, z = 1), evaluated for p = p m>k = 
^|0m,fe| 2 , and Ifc is introduced accordingly. Expressing 
the O(N) symmetry in the 2-replica second-order tensors 
is a little more tedious, but nonetheless straighforward. 

The resulting flow equations in scaled form read (where 
for ease of notation we drop the subscript k in the right- 
hand sides, i.e., up to a sign, the beta functions, for 
all quantities but T k and also drop the argument of 

v{pi,P2,z))- 

d t Uk(p) = -(d-2 + fj- Jj)u(p) + {d-4 + rj) pu'(p) 
+ 2vJ(N - l)l{ d \w T (p))S T (p) + l[ d \w L {p))8 L {p) 



+ 2T k v d \ (N - l)i d \w T (p)) + l { d \w L (p)) 



(77) 



d t v k (pi,p 2 , z) = -(d-A + 2fj- 2rj)v + (d-i + fj) 



(pid pl +p 2 d P2 )v--^{(N-l) 



ip 2 ii(w T (pi)) 



5 T (p 1 )(2p 1 d Pl v - zd z v) + l[ d l(w T (p 1 ),w T (p 2 ))(d z v) 
+ (1 - z 2 ) [4p 2 l 2 d \w T (pi))S T (pi)d 2 z v+ 
8pll[ d l{w T {pi),w L (p 2 ))(d P2 d z v) 2 - l{ d l(w T (pi),'w T {p 2 )) 
((d z v) 2 + 2zd z vd 2 v - (1 - z 2 )(d 2 v) 2 ) 
+ 8 PlP2 \l {d \w L ( Pl ))S L (p 1 )(d Pl v + 2 Pl d 2 p v) 



2 PlP2 l (d l{w L {pi),w L { P2 )){d Pl d P2 v) 2 

-T k ^{(N- l)p2lf(w T (p 1 ))(2p 1 d Pl v - zd z v)+ 
P1P2 I- 

(1 - Z 2 ) P 2li(w T ( Pl ))d 2 z V + 2 PlP 2l d 1 (w L (p 1 )) 

(d pi v + 2pid 2 pi v)+perm(l2)}, 



perm 



(12)} 



(78) 



Vd ) V 



m ( 2 l{w L {p m ),Q) +m { 3 d l(w L {p m ),Q) 
w L {p m ) 2 



$T(p m ) 



+ 8m (d l(w L (p m ),0) (79) 



Pn 



x w L (p m ) 



SriPm) - S L (p m ) 
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2l]k - Vk 



2v d 



PmU"(p m ) 

ll(0)u"( Pm ) + l 2 d \w L (p m ))S L (p m ) 



(N - l)p m lf(0)6' T (p m ) 2 + 

(1 + 2p m 5' T (fi m ) 



+ 2p 2 m 6'+(p m ))u"(p m ) - 2p 2 n u"'(p m )5' T (p m ) 



21^1(0, WL{Pm))u" \p m ) 



1 + PmS' T (p r , 



p m lf ] ' {w L (p m ))5' T (p m )5' L (p m ) \ + 



(80) 



where all symbols have the same meaning as in the 
previous equations and, by construction, Wi(p m ) = 
2p m u"(p m ), WT{p m ) = 0, and Sripm) = 1- Note that 
in the last two equations, we have omitted for simplicity 
the (subdominant) terms involving T k in the beta func- 
tions and that in Eq. (|80p. the dots denote "anomalous" 
terms which generalize those found for the RFIM (see 
Eq. (|69p ) and vanish when the function Vk(pi, Pi, z) is 
analytic in all its arguments; their expression is lengthy 
and will be discussed in the companion paper^ 

When N = 1 and z = ±1, Eqs. tfTTJ and (JTSJ) re- 
duce to the previous equations for the RFIM, Eqs. 
and (|62j). expressed with p as variable instead of 0: 
Vk(pi,P2, z = +1) is equal to v k ((fii, ^2) for ^1^2 > and 
v k (pi,P2,z = -1) is equal to v k ((pi,<p2) for ipiip 2 < 0; 
5 k ,h{p) = S k ((p) and w k x{p) = u"((p)M Finally, the 
comments made about the important features of the flow 
equations for the RFIM carry over to the equations for 
the RFO(7V)M. 



D. Application to related disordered models 

Even though we have chosen to more specifically focus 
on the random field model, it is worth sketching at this 
point the relevance of the NP-FRG equations derived in 
this section to other disordered systems. (As stressed 
already several times, we exclude spin glass ordering from 
our considerations.) 

The flow equations obtained for the RFO(A)M, 
Eqs. ((7711751), directly apply to the RAO(A)M for de- 
scribing the long-distance physics associated with ferro- 
magnetic ordering. The putative fixed points are also 
expected to be at zero temperature, so that similar scal- 
ing dimensions need be introduced. The specificity of the 
random anisotropy model comes in the initial conditions 
(see section II-A) and in the additional symmetry of the 
2-replica potential, namely, Vk{pi,P2,z) = Vk{p\, P2, — z). 

Similarly, the flow equations for the RFIM, 
Eqs. (]61l62l66p . can be applied to the random elastic 
model. In this case, one can check that, owing to the 
statistical tilt symmetry, u' k (tp) = and r\ k = while 
v k (tpi, ^2) = Vk{<fi — ^2) • After introducing the variable 
y = (pi — ip 2 and dropping the temperature, Eq. (]62p . 



can be rewritten as 



-d t v k {y) = (d-A + 2fj k )v k (y) - -(d-4 + fj k )yv' k (y) 



2v/ 2 d \0) K(y)-2<(0)K'(y), 



(81) 



where a prime denotes a derivative with respect to y. The 
roughness exponent is defined through ( = — (d— 4+?y)/2, 
and one can then see that the above equation reduces 
to the one-loop FRG equation for a disordered elas- 
tic medium i 28 ! 44 Going beyond this level of description 
requires to consider the next orders of the truncation 
scheme, in particular to include the 3-replica potential 
and apply the next order of the derivative expansion for 
the 2-replica effective average action. 

Finally, Eqs. (|61I62I66[) can be used in the case of the 
random temperature model with an appropriate account 
of the symmetry: u k = u k {p), v k = v k (pi,p2), with 
p = (p 2 /2. However, the scaling dimensions introduced to 
search for a zero-temperature fixed point are not appro- 
priate in the present case where one anticipates a fixed 
point at a nonzero temperature (for a preliminary non- 
perturbative treatment, see Ref. [68|). 



RECOVERING THE PERTURBATIVE 
RESULTS 



Analysis of the NP-FRG equations near d = 6 
and for N — > 00 



For ease of notation, we only consider the RFIM, but 
a similar analysis holds for the RFO(7V)M. It is easy to 
check that the flow equations, Eqs. (pTl I62I66I69|) . admit 
for fixed-point solution the Gaussian fixed point charac- 
terized by V i G) = fji G) = 0, u[ G) (<p) = 2v d l[ d) (0)/(d- 2), 

and (5; (1^1,(^2) = 1- The Gaussian fixed point is once 
unstable for dimensions larger than 6, but the coupling 
constant associated with the c/2 4 -term in u((p) also be- 
comes relevant for dimensions less than 6 so that the 
Gaussian fixed point becomes unstable for d < 6, as al- 
ready well known. 

The first order in e = 6 — d can be derived by a di- 
rect expansion of the fixed-point solution, with u*(tp) = 

(G) 

ul +eu 1 (ip), 5*((pi,<p2) = 1 + eSx(<pi, ¥2)- One easily 
finds that at this order one still has 77* = fy» = 0. After 
inserting these results in Eqs. (|61I62[) . deriving the equa- 
tion for Vk with respect to (pi and ip 2 , and setting the 
left-hand sides to zero, one obtains the following equa- 
tions for u\(ip) and S\(ipi, (p 2 ): 



= 4 Ul (<p) 



11 \ Ve 



4 6) (0)-44 6) (0KM 



+ 4^ 6) (0)<5 1 (^) 



(82) 
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= {ip 1 d v>t +ip 2 d lp2 )8 1 {ip 1 ,^2) 



2v 6 l 2 6) (0)(d Vl +d V2 ) 2 6 1 (p ll p 2 ). 



(83) 



Vm*) 2 + constant with 



By introducing the variables x = (pi + p2)/2 and y — 
(pi — <f2)/2, the latter equation can be rewritten as 

(xd x + yd y )S 1 (x,y) = 2v 6 l 2 6) (0)dlS 1 (x 1 y). (84) 

The symmetry of v(pi, ip 2 ) with respect to the exchange 
of ipi and f2 and to changes of sign of ip\ and if 2 (see IV- 
B) translates into the fact that Si is an even function of 
x and y. Provided one requires that Si(x, y — » 0) is finite 
(which is needed for a well defined renormalizable the- 
ory) , the only acceptable solution satisfying this property 
is a constant; due to the constraint 5 m ,k = Sk(p>m,k) = 1, 
it is equal to zero, i.e., 5%(x,y) = 0. In addition to the 
now unstable Gaussian fixed point, Eq. (1521 has then for 
solution ui(ip) = (Xi*/8)(p 2 

(Pm* = 6v e l 2 6 \o). One also finds Ai* 
so that, up to irrelevant constant factors, the solution 
corresponds to the fixed point of the pure Ising model 
(no random field) at first order in e = 4 — d. The fixed 
point is found once unstable and the associated expo- 
nents, e.g., v = i + satisfy the d — > d— 2 dimensional 
reduction. 

Equivalently, one can make a more direct connection 
to standard perturbation analysis by reframing the above 
results in a double expansion in e and in the p A coupling 
constant defined through Afe = u' k ' (p m ^ k ). Introducing 
as before p m ^ k = (1/2)^ k , one obtains from Eqs. ([771 - 
[50]) th at 77,77 = 0(X 2 ), 5 = l + 0(X 2 ) and 

d t \ k = -eX k + 36v d l {d) (0)X 2 k + 0(XlT k Xl), (85) 



dtPr, 



-(2-e)p m}k + 6v d [l 2 d) {0) 



44 d) (0)A fePm , fe ] 



0(A 3 fc ,T fc A 2 fc ), 



(86) 



where we have used the Taylor expansion of the threshold 
functions for small arguments. (The fixed-point solution 
of Eqs. (|85l86p is of course equal to that obtained above 
with A* = eAi* and p mit = p 2 ntf /2.) Again, up to irrele- 
vant factors, this gives back the one-loop perturbative re- 
sult for the pure Ising model obtained in a weak-coupling 
expansion in d = 4 — e. 

The above result is derived through an expansion in a 
single coupling constant, Afe, associated to the 1-replica 
part of the effective action. It has been argued by Brezin 
and De Dominici a 69 ' 70 that one should consider instead 
an expansion involving all ip 4 coupling constants associ- 
ated with multiple replicas. In the present formalism, we 
can perform a more careful analysis using the tp 4 cou- 
pling constants associated with the 2-replica part of the 
effective action, coupling constants that are considered 
as potentially relevant in Refs. [69.70]. We find that this 
does not change the conclusion and, as previously ob- 
tained in Ref. [7l| , that the fixed point corresponding to 



dimensional reduction is still once unstable at first order 
in e. This is discussed in more detail in Appendix A. 

The above analysis is extended to the O(N) version in 
a straightforward way. The property that the perturba- 
tive result at first order in e = 6 — d is recovered within 
our nonperturbative approximation scheme is actually a 
consequence of the one-loop-like structure of the exact 
flow equation for the effective average action, Eq. (|16p . 
For the very same reason, the large N limit can also be 
easily recovered. 

Rescaling the variables as p — » Np, z — » z and the po- 
tentials as u — * Nu, v — > Nv, and retaining only the dom- 
inant terms when N — > 00, one finds that 77 = 0(1 /N), 
fj = 0(1 /N) and that the "longitudinal" contributions 
drop out from the RG flow equations. As a consequence, 
Eqs. ([77)) and ([75]) can be recast as 



d t u k (p) = -(d - 2)u k (p) + (d- 4)pu' k (p) 



2v d l[ d \u' k (p))S k , T (p), 



(87) 



d t Sk,T(pi,P2,z) = (d-<i)(pid pi + p2d P2 )5 k ,T(pi, P2, z)- 

l[f«(pi),U> h (p2)) 



2v d 



y/Plp2 

4V*(pi)) 



-Sk,r(pi,P2, z)d z 5k,T(pi,P2, z) 



2 P i 



h,T(pi)(2pid Pl - zd z )S k ,T(fii, P2, z) 



2 P 2 



$k,T(p2)(2p2d P2 - zd z )8k,T(pi,p2,z) >, 



where we have defined a generalized "transverse" disor- 
der cumulant S kj x(pi, P2, z) via an extension of Eq. ([74p . 
namely, 



$k,T(pi,P2,z) 



2^P\P2 



d z Vk(pi,P2,z), (89) 



which reduces to 5 k ^T(p) when p\ = p 2 = p and z = 1. 
Eq. ([55)) is obtained by deriving the flow equation for 
Vk(pi, p2,z). 

If one starts the flow equations with an initial condition 
va(pi,P2,z) = 2y/pxp 2 z (corresponding to S a ,t = 1), 
the beta function is identically zero and one therefore 
finds that the solution of Eq. ([88]) at all scales remains 
Sk,r(pi, P2, z) = The resulting equation for the 1- 
replica potential is then very similar to its counterpart for 
the pure O(N) model with N — > 00 limit in dimension d— 
2 (the flow equation is then simply given by the LPA^) . 

To see more explicitly the connection, one can follow 
the flow of the (p 4 coupling constant Afe = u'l(p m _ k ) as 
well as that of p m ,k which, we recall, satisfies u' k (p m _ k ) — 
and is akin to a (dimensionless) order parameter at the 
running scale k. One finds 



-(6-d)A fc 4 
+ [(d 



4v4 d) (0)Xl 



2v d l {d) ( 



'(Pm,k), 



(90) 
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d t p v 



-(d-4)p TO)fc + 2^4 d) (0) J 



(91) Eqs. d79j) gives 



which results in the nontrivial fixed point p m * = 
2v d l ( 2 d) (0)/(d - 4), A* = (6 - d)/(Av d l [d) (0)). This fixed 
point is once unstable (and it remains so when consider- 
ing the additional directions associated with the 2-replica 
potential, see above) and is characterized by critical ex- 
ponents satisfying the dimension reduction property, e.g., 
v = l/(d — 4) to be compared to v = l/(d — 2) for the 
pure model. 

Note that the above perturbative expressions are re- 
covered from the truncated NP-FRG equations even 
with an additional approximation using a field expansion 
around the minimum of the 1-replica potential. 



B. Recovering the perturbative FRG near D = 4 

A strong property of the minimal nonperturbative 
truncation described above is that it also reduces, in the 
appropriate limit and for the RFO(iV > 1)M, to the per- 
turbative FRG equations at first order in e = d— 4 derived 
by Fisher £2- The latter are obtained from a low-disorder 
loop expansion of the nonlinear sigma model associated 
with the RFO(iV)M. It is therefore quite remarkable that 
our formalism in which no hard constraint is enforced 
leads to the proper result within the minimal approxi- 
mation scheme. 

For the RFO(7V)M with N > 1, d = 4 is the lower 
critical dimension for ferromagnetism. (We mean here 
long-range ferromagnetic order with a nonzero order pa- 
rameter, the case of quasi-long range order will be dis- 
cussed later on.) As a result, the critical point and the 
associated fixed point occur near d = 4 for a value of p m 
that diverges as 1/e with e = d — 4. As in the case of 
the pure O(N) model near d = 2)21 one can therefore 
organize a systematic expansion in powers of l/p m - 

At the minimum of the 1-replica potential (p = 
p m ), the transverse mass, associated with the Goldstone 
modes, is zero whereas the longitudinal mass is very large 
and scales as p m (anticipating that u"(p m ) does not van- 
ish). One can then use the asymptotic properties of the 
threshold functions for large arguments, 



lid) 



(w 



oo 



w 



-(m+l) 



m 



(d) 



, (w — ► oo, 0) 



w 



(92) 
(93) 



which encodes the decoupling of the massive mode. 

In addition, we assume that as p m — ► oo, o~L,T(Pm) 
stay finite (recall that actually, SriPm) — 1) and that 
their derivatives, S' L T (p m ), etc, go to zero at least as fast 
as l/pm] on the other hand, p m is a singular point for 
u(p) (the location of its minimum), so that even when 
we expect that u"(p m ), u"'(p m ), etc, stay of 0(1). The 
consistency of these assumptions is easily checked a pos- 
teriori. Inserting the above results and assumptions in 



8v d 



dp. 



m.k 



(94) 



which shows that r\ is of order 1/ p m . 

Deriving once the flow equation for the 1-replica po- 
tential Uk(p) leads to 

dtu' k (p)=-(2- Vk )u' k (p) + (e + fj k )pu , /(p) 

-2v d (N-l)i d \u' k (p))5 KT (p), (95) 

from which one obtains the flow equation for the running 
order parameter p m ,k- 



dtpr, 



-(e + fj k )p m , k + 2v d (N - l)4 d) (0), (96) 



where e = d — 4. (Note that we have again omitted the 
subscript k in the right-hand sides and dropped the sub- 
dominant terms involving the renormalized temperature 
T k .) The last equation shows that the fixed point value 
of p m ,k satisfies, as anticipated, p m * = 0(1/ e), which 
results in 77, 77 = 0(e). 

One can now apply a similar treatment to the flow 
equation for the 2-replica potential evaluated for p\ — 
P2 = Pmk- For convenience, we introduce the function 



Rk(z) 



z ) 



(2p 



m,k) 



(97) 



which, due to Eq. (|74l) and the constraint Sk,T(pm,k) = 1, 
satisfies R' k (z = 1) = l/(2p TOi fc).— The flow equation for 
Rk(z) can be expressed as 



d t R k (z) 



— -^dtv k (p,p,z)\ p=Pn 

\^Pm,k J 



dtPm,k dp 



Vk(p,P,z) 



(2p) 2 



\P=P-n 



(98) 



which with the help of Eq. ([96)) finally leads to 



d t R k (z) ~ (e + 2 Vk )R(z) - 2v d li d \0){(N - 1) 
AR k (z)R' k (l) + R' k (z)(R' k (z) - 2zR' k (l)) 



(1 



R' k (z) 2 + 2(R' k (l) - zR' k (z))R'/(z) 



+ (l-z 2 )R'/(zf 



(99) 



To dominant order in e, one can set d = 4 in v d and 
4^(0) m au equations and in v d /d in Eq. Using the 
property of the threshold function 1% (0) = 1 + 0(rj) 
and discarding subdominant terms, one finally arrives at 

Vk = iv 4 R' k (l), f} k = -e + 4(N - l)v 4 R' k (l), (100) 
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d t R k (z) =eR k (z) - 2v 4 ^A(N - 2)R' k (l)R k (z) + 



(N-l) 



R' k (z)-2zR' k (l)) 



R' k {z) 



+ (1 - - 2 ) 



(101) 



-R' k (zf + 2(R' k {l)- 



zR' k {z))R' k \z) + (l-z 2 )R' k \zf 



where 



= 32-7T 2 and R k (z) is of order e near its 



fixed point. The above equations coincide with the one- 
loop perturbative FRG equations derived by Fisher^ 
Note that this result is independent of the choice of 
the infrared cut-off function R k (q 2 ): indeed, one eas- 
ily checks that not only 1^(0) = 1 + 0(rj), but also 
lim^^oo [m^liw, 0)u; 2 ] = 1 + 0(rf), irrespective of the 
regulator. 

Finally, we note that setting N = 2 and introducing 
the variable <f> — cos^ 1 (z) in Eq. (jlOip leads to 

- d t R h {<j>) = eR k (4>) - 2v 4 [R'l{<j>) - 2J#(0)] R'l{<j>), 

(102) 

which, after use of Eq. p00[) for r] k and rj kl coincides 
with the 1-loop perturbative FRG equation for a disor- 
dered periodic elastic system with a one-component dis- 
placement field: compare for instance with Eq. (|81[) . in 
which one should set £ = due to the periodicity. 72 (Be 
careful, however, that rj k and fj k denote different sets of 
exponents in the formalism leading to Eq. ([8T]) and in 
the present one.) 7 ^ 



VI. CONCLUDING REMARKS 

In this work, described in the present paper and in the 
following one^Si we have developed a theoretical approach 
which is able to describe the long-distance physics, crit- 
ically, phase ordering or "quasi" -ordering, of systems in 
the presence of quenched disorder, in particular random 
field models for which standard perturbation theory is 
known to fail. 

The approach is based on an exact renormalization 
group equation for the effective average action (the gen- 
erating functional of 1 — PI vertices) and on a nonper- 
turbative truncation scheme. This nonperturbative RG 
formalism has recently been applied with success to a va- 
riety of systems. The key point in the present problem 
is to provide a proper account of the renormalized dis- 
tribution of the quenched disorder, and we have shown 
that this can be conveniently done through a cumulant 
expansion and the use of a replica method in which the 
permutational symmetry among replicas is explicitly bro- 
ken. 

We have stressed that any relevant treatment of ran- 
dom field models and related disordered systems must 
include the second cumulant of the renormalized disor- 
der, i.e., at least a function of two (replica) field argu- 



ments. Accordingly, we have proposed a nonperturbative 
approximation scheme. Within this scheme, the min- 
imal truncation for the RFO(iV)M already reproduces 
the leading results of perturbative RG analyses near the 
upper critical dimension, d uc = 6 and when the number 
of components N becomes infinite. More importantly, 
it gives back the perturbative FRG equations near the 
lower critical dimension for ferromagnetism when AT > 1 , 
d = 4. 

One of the main advantages of the present approach, 
which will be illustrated in the following paper, is that it 
provides a unified framework to describe models in any 
spatial dimension d and for any number N of field com- 
ponents. As such, it garantees a consistent interpolation 
of all known results in the whole (N, d) plane, in ad- 
dition to allowing the study of genuine nonperturbative 
phenomena. 

We thank D. Mouhanna for helpful discussions. 



APPENDIX A: EXPANSION IN SEVERAL 
COUPLING CONSTANTS NEAR d = 6 

Near the upper critical dimension d uc = 6, the flow 
equations for the RFIM derived within the minimal non- 
perturbative truncation, Eqs. (|6T1 l62l [66l 169]) , can be 
expanded in several cp* coupling constants, in order to 
make the connection with recent one-loop studies of the 
RFI M 69 i 70 i 71 . On top of the 1-replica ip 4 coupling con- 
stant already used in section V-A, X k = u'l" (ip m ,k) , we in- 
troduce two additional coupling constants obtained from 
the 2-replica potential, 



1 



u k,2 = -- {d 2 Vl + d* 2 ) d Vl d V2 v k {ip 1 ,p 2 )\ ipi = v . 1 = 



(Al) 

u k ,3 = -dl 1 dl 2 v k (tf 1 ,Lp 2 )\ Vl = V2 = lfmkl (A2) 

which amounts to consider a 2-replica potential of the 
form 



Uk,3 



U k ,2 



{<Pl<P2 ~ 4<Pm,k) + 



(A3) 



where the dots denote higher-order terms in the field ex- 
pansion around the minimum of the 1-replica potential 
and, as before, S„ hk = d ipi d ip2 v k (tpi,ip 2 )\ Vm k = 1. The 
present description is thus very similar to that used in 
Refs. [69:.70|, except that we do not include 3- and 4- 
replica terms. However, the issues raised by Brezin and 
De Dominici o 69 ' 70 can already be addressed by consider- 
ing the 2-replica term. 

By expanding the flow equations for the 1- and 2- 
replica potentials in powers of the coupling constants, 
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which we generically denote u kia with u k \ — X k , one 
finds that r\ = 0(u 2 ) and that up to a 0(u 3 ), 

dtu k s = (d- 6)u M + 2-y d |6Z^ ) (0)^, :1 + 124 d) (0)w fc ,iX 

(A4) 

9 t Wfe,2 = (d- 4)M fej2 + 2v d \ 6l 3 d \o)u kA (u k: 2 + U k ,3)+ 



with, we recall, p m , k — tp 2 m k /2. In addition, by using 
d t T k = (2 + % - »7k)Tfc and the equation for 2r) k - fj k , 
Eq. ([69]) . one obtains to a 0(m 2 ): 



d t T k = 2T k + Av d T k { 2l ( 2 d) {0){u k:2 + u k , 3 ) + 6l[ d) {0)x 



+ T k l\ '{Q)(2u k ^ + 3u k , 3 ) 



U k A 



(A8) 



6l 2 d \o)u k:2 (u kt 2 + u kt3 ) + 3T k l 2 a >(0)u kA u k: < 

(A5) 

d t u k>3 = (d -4)u fe ,3 + 2t> d |4 <0 (0)u^ il +44 d) (0K,ix 
(u k ,2 + u kj3 ) + 2l 2 d \0) [(uk,2 + u kt3 ) 2 + 3u 2 k3 ] 



(d), 



+ 2T k l < 2 d) {Q)u k>1 u k<3 \, 



(A6) 



A/c/3 m ,fe 



2Z^(0)K, 2 + u M ) + T k l\ ai (0)u kjl + 



24 d) (o)u fc ,i + 8i ( 2 a) (o)(u k , 2 + «M) + ^ flJ (QK 



Afe/3 m ,fe 



(««)/ 



It can now be checked that the above equations co- 
incide with those derived by Mukaida and Sakamoto 71 
with the introduction of new running coupling constants: 
.go = Tfe(5fe(0) _1 ,5i = u k ,iSk(0),g2 = u k ^,g 3 = u k , 3 , with 
4(0) = S k (cpi = 0, (p 2 = 0) = l + 2(u ki2 + u k ^)p mtk (and, 
of course, the 3- and 4-replica contributions missing). As 
these authors, we therefore obtain that the dimensional 
reduction fixed point is once unstable at first order in 
e = 6 - d. 

On the other hand, the analysis performed by Brezin 
and Dc Dominici a 69 i 70 requires to introduce different scal- 
ing dimensions, corresponding to new coupling constants, 

9o = 9o,9i = 51)52,3 = 92,39a 1 - The beta functions 
we now obtain for go , g\ and 92 coincide with those of 
Ref. [69J7OJ, but that for g 3 is ill defined as it contains 
a term that blows up as k — > 0. The scaling dimensions 
suggested by Brezin and De Dominicis are thus not com- 
patible with our approach. 
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